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Abstract 

Binary nonlinearization of AKNS spectral problem is extended to the cases of higher-order symmetry 
constraints. The Hamiltonian structures, Lax representations, r-matrices and integrals of motion 
in involution are explicitly proposed for the resulting constrained systems in the cases of the first 
four orders. The obtained integrals of motion are proved to be functionally independent and thus 
s ! ' the constrained systems are completely integrable in the Liouville sense. 



O 1 Introduction 



Symmetry constraints have been applied to constructing finite-dimensional integrable systems from Lax pairs 
of soliton equations. It has aroused an increasing interest among scientific researchers in the field of soliton 
theory in the past several years. A large class of finite-dimensional integrable systems had been obtained 
indeed from soliton hierarchies by imposing symmetry constraints 0~lH. There is also a fc-constraint KP the- 
ory which uses symmetry constraints to extend the standard KP theory Nonlinearization methodo 
utilizes a kind of special symmetry constraints for 1 + 1 dimensional soliton equations. In refs. 

0-0 

, one 

of the authors suggested binary nonlinearization method to get finite-dimensional integrable system from 
Lax pairs of 1+1 dimensional integrable systems. This method has been successfully applied to a few of fa- 
mous soliton hierarchies_and it has to be adopted in making nonlinearization when Lax pairs have odd-order 
matrix representations^ . 

However, what has been usually considered in binary nonlinearization is the Bargmann symmetry con- 
straints (explicit constraints). In this paper, we would like to extend binary nonlinearization method to 
the cases of higher-order symmetry constraints (implicit constraints). The paper is organized as follows. In 
Section 2, the Hamiltonian structures are constructed for the resulting constrained systems. Main analysis 
is to describe how to introduce new dependent variables to put the constrained systems to be Hamiltonian. 
In Section 3, the same r-matrices as those in mono- nonlinearization are proved to be valid for the con- 
strained systems in binary nonlinearization after presenting Lax pairs. In Section 4, the obtained integrals 
of motion are shown to be in involution in pairs and functionally independent. It follows that all those 
finite-dimensional dynamical systems are Liouville integrable. A conclusion and some remarks are given in 
the final section. 

2 The Hamiltonian structures 

For the sake of simplicity, we assume that the potentials v, w and all their derivatives with respect to x tend 
to zero when \x\ — > oo. 
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Following the notation in the referencc01 , we consider the spectral problem 

{ 1 x i.e. <t> x = U<j>, = \,U=U(u,\)=[ , ,« = (2- 1 ) 

[ 4>2x = W(pi + A0 2 , V^ 2 / \ W * J \ W ) 

and its adjoint spectral problem 

f 1pi X = Xlpl - W1p 2 , rp ( \ 

i.e. V, = ^ = , (2.2) 

[ i>2x = -Vlpl ~ \lp2i \ W2 ) 

for the AKNS soliton hierarchy. Assume that V x = [U, V] has a formal solution 

v=[ a b ) = y(<* bi )A _ 
c -a I \ Ci -en 

' i=Q v 

where cij, &i, Cj, z > 0, are determined uniquely by 

= qci - r6j, & ix = ~2b i+ i - 2gOj, c ix = 2c i+ i + 2raj, i > 0, 

and 

a = -1, b Q = c = 0; a 4 |„ =0 = 6j|«=o = c 4 | M=0 = 0, i > 1. 
Upon making the Bargmann constraint, as in ref. 0, 



h \ _ ( v \ _ ( <i) 2 ,(f>i> \ _ ( J2iLi ^2 Ai)0i(x,Ai 



(2.3) 



where we use 0j, ipj, j = 1,2, to denote the following vector valued functions 

<t>j = (<Aj(£ I Ai),...,0 i (x,AAr)) T , Vj = (^j(x,Xi),-;ip](x,XN)) T , 3 = 1,2, 
we obtain from (2.1) and (2.2) the following Hamiltonian system 

P =-°* Q (24) 

x oq ' dp 1 ' 

with the Hamiltonian 

i? =< Aip2,<h > < Atpi,4>i > + < -02,01 >< -01,02 >, v4 = diag(\ u ...,\ N ), Aj ^ A^ (i ^ j), (2.5) 
where 

Q = (0i (x, Ai), ...,0i(x, Aat),0 2 (x, Ai), 2 (x, Xn)) T , 
P = (ipi(x, Ai), ...,ipi(x, Ajv), ^(x, Ai), -02 (x, Aat)) t . 

We will focus on the higher-order symmetry constraint 

c n \ SHn-l /<-0i,0 2 >\ 2 

&„ J := -*r =/ H <^ 2 ,0i> J' Hn - i= n a - +1 ' (2 - 6) 

where /3 is a constant to be determined and n is a positive integer. If n is odd, we can introduce the 
Jacobi-Ostrogradsky coordinates to put the constrained system consisting of (2.1) and (2.2) with all A = Aj 
and (2.6) to be Hamiltonian. For example, in the case of n = 3 we can introduce new dependent variables 

0JV+1 = V, 0W+2 = W, -0AT+1 = --0AT+2,i, -0JV+2 = ~^4>N+l,x- (2.7) 
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The constraint (2.6) with /3q = 1 becomes 



C3 



1 / W xx - 2VW 2 \ _ ( - ( l ) N+2,x - \<j)N+l4> 2 N+2 \ / < ^1,02 > 



h ) A\ v xx -2v 2 w J I , 1j2 j I V < -02,01 > 

1 — <?N+2,x - tjVn+iVN+Z ' 



(2.8) 



We take the Hamiltonian H for the constrained system as follows 

H =< Ai> 2l 4>2 > - < Alf)x,(j)l > +<j>N+l < 01,02 > +4>N+2 < "02,01 > + ^ 

and thus the constrained system can be written as 



4>li,x = -g^- = — Ai^li + 0AT+102i, 4>2i,x = -gn^: = -402i + 0AT+201J, 

OH 3H 

= 30^1 = ~Wn+2, 0N+2,x = g^r 2 = -4lpN+l, 

f)ff (9/7 

^JV+l.x = - ^ N+1 = - < "01,02 > -20!r +1 0jV+2, VaT+^s = - ^ jy+2 = - < "02,01 > ~ 20^ 0W+2 , 

(2.9) 

where 0^ = (f>j(x,\i), ipji = ipj(x,\i), j = 1, 2, 1 < i < N. 

In general, when n = 2fc + 1, the symmetry constraint (2.6) with /3 = 1 is 



C2fc+1 \ _ 8H 2 k ( < -01,02 > 



&2fc+l / 5u \ < -02, 01 > 
Introduce the following Jacobi-Ostrogradsky coordinates!! 

i d^V i d l ~ lr W ' -i n J, 

07V-l+2^ — ^3=1 i 0JV+2i = i I = 1, 2, fc, 

where 

dx l dx 1 ' 8v® f^ ' 0uC»+O ' ^ v > dw^ dx' - 

and denote the dependent variables by 



(2.10) 



Q — (011, 017V, 021, 02JV, 0AT+1) 0W+2fe) T , 
P = (-0H, 0^,-021, ...,02AT,0W+1, ...,0AT+2fe) T - 

The constrained system consisting of (2.1) and (2.2) with all A = A^ and (2.10) becomes 



(2.11) 



= _dH2k ±1 = dH2k+i . . 

dQ ' ^ dP y ' 



with the Hamiltonian 

2fe 



H2k+l — E 4>N+i,xl/->N+i — #2fc + H 2 k+l, 



1 = 1 
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where the third term is given by 

H 2k +l =< Alp 2l (f>2 > - < Atpi, 01 > + 0JV+1 < "01, 02 > + 0W+2 < "02, 01 > • ( 2 -13) 

When 0ii, tpn, 02i, tp2i (i = 1, —,N) satisfy (2.1) and (2.2) with A = A,, we have (see (3.18) in ref. 0) 
k =< A l - 2k -^ 2 , 0! >, q =< A l ~ 2k - l i> u 2 >, ai = i(< A'- 2fc - Vi, 0i > - < A 1 " 2 *- 1 ^, 02 >), (2.14) 

which are used in generating integrals of motion. Therefore all the constrained systems possess the Hamil- 
tonian structures. 

When n is even, the regular condition for introducing the Jacobi-Ostrogradsky coordinates can not be 
satisfied . Nevertheless we still can introduce new dependent variables but not the Jacobi-Ostrogradsky 
coordinates, to put the constrained system to be Hamiltonian. For example, in the case of n — 2, since 
b 2 = —T£V x ,c 2 = the constraint (2.6) with 0q — — ^ reads as 

V x =< 1p2, 01 > W x = - < Ipl, 02 > ■ (2-15) 

Now we do the trick by introducing new dependent variables 

0AT + i = V, 4>N+1 = W. 

Two vectors of all dependent variables are given by 

Q = (011, 01JV, 021, 02W, 0Af+l) T , P = ("011, ■■;i>lN,1p2l, ■ ■ ■ , "02JV , "0JV+1 ) T ■ 

It is easy to verify that the constrained system can be transformed into the following Hamiltonian form 



F) T-f FiT-f FiT-f 

hi,x = -g^p = —Xilpli + 0AT+102i, <j> 2 i,x = -g^- = Ai02i + 0JV+101i, 4>N + l,x = Q^ N+1 = < ^2, 01 

FiT-f f)TJ f)TJ 

hi,x = = ^iV'li - 0AT+l'02i, 1p2i,x = = ~ 4>N+l^lU tpN+l,x = = ~ < V 

(2.16) 

where 1 < i < N and the Hamiltonian is 



H =< A0 2 ,02 > - < A0l,0l > +0AT+1 < "01,02 > +'0JV+1 < ll>2, 01 > • 

In the case of n = 4, since 6 4 = ^(—v xxx + 6vv x w), C4 = ^(w xxx — 6vww x ), the constraint (2.6) with 
flo = — \ reads as 

\{v xxx - 6vwv x ) =< ^2,01 >, j(^xxx - 6vww x ) = - < ipufa > ■ (2-17) 
Introducing new dependent variables 



1 IS 

07V+1 = V, 0JV+2 = —jVx, 0AT+3 = sV xx - %V 2 W, 

1pN+l = \w xx - |l>U> 2 , ^JV+2 = \w x , TpN+3 = W, 



(2.18) 



and two vectors of all dependent variables 

Q = (011, 01JV, 021, 02N, 0JV+1, 0W+2, 07V+3) T , P = C011 , ■ ■■ , IplN , "021 , ■■ ; lp2N , IpN+l , 1pN+2 , i'N+zf , 
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the corresponding constrained system becomes 



<f>li,x = -AjVli + 4>N+l<t>2i, 
<j>2i,x — \4>2i + (f>N+l<Pli, 
4>N+\,x = —%4>N+2, 

<t>N+2,x = - ^(j) 2 N+1 1p N+3 , 

4>N+3,x =< lfa,4>l > -|^+i^JV+2 - \<t>N+l<t>N+2^N+3, 
tpli,x = AiV'H - 4>N+3tp2i, 
1p2i,x = -\lp2i - 4>N+l1pli, 

IpN+l.x — ~ < 1pl,<f>2 > + ^<l>N+l<j>N+2<pN+3 + §0^+2^+3 
1pN+2.x = f^iV+l^ir+3 + TpN+1, 



(2.19) 



k "0^+3,2; = 2lp N+2 , 



where 1 < i < N. We can find that the Hamiltonian for this system (2.19), which can be selected as follows 



We note that when n = 2 or 4, the new dependent variables are some kind of combinations of b\,C\ 
or bi,b 2 ,b 3 ;ci,c 2 ,c 3 . In general, for the case of n = 2k, new dependent variables may similarly be taken 
as combinations of b\, &2fc-i! ci, c 2 k-i- However we do not have any concrete expression for such new 
dependent variables. How to introduce them still needs further investigation. 

3 Lax representations and r-matrices 

For the system 



H = < Alp 2 , <t>2 > ~ < > +<t>N+l <1pl,<h> +1pN+l < 1p2, 01 > 

-20N+20AT+1 - 2<f> N+3 1p N+2 - ^(f> 2 N+l 1pN+2lpN+3 ~ | 0JV+10JV+2^jv+3 ■ 



(2.20) 





(3.1) 



we take the symmetry constraint 




(3.2) 



where the /3q is a special constant for each case of n. Assume that 




(3.3) 



where 



2 (^lj <t>lj ~ 4>2j <f>2j ) , bj= 1p2j <f>lj , Cj = 



i < i <n. 



At this moment, we have 




(3.4) 
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The reasons are the following. We first have 

N f ~ 

N x n) = \4 (n) + ^0* = V^ n) + ( °J X 

3 = 1 

N 



A - A ., 



[u, *w ] = [c/, v w ] + A,£ w C J ~ vb l ~ 2 b ^ - 2v i a > 

; , A A • V - r " .. +2A Cj -w Cj +v bj 

N 1 / ~ ~ ~ ~ \ 

= [U, F(«)] + A £ T 2 !- ^ ~" ^ " 2 6j ^ " 2 " ° J + A 



w Cj — v bj —2 foj Xj — 2w a,j 
Then taking the residue of (3.4) at A = Xj , we get the equations 



-2 < V 2 ,0i > 

2<^ 1 ,^a> 



4>ij{i>2jx + AjV'aj + W'y) + ^2j(4>ijx + ^]4>ij - w^i) = o, 

4>2](^ljx ~ Xjlplj + W1p2j) + 1plj(4>2jx ~ Aj^2j ~ W'lj) = 0, 
~4>2j{lp2jx + Xjlfaj + Vlp 2 j) + 1plj(<t>ljx + Aj^lj - W^j) = 0, 
4>lj{^ljx - Xji>lj + Wlp2 3 ) ~ 1p2j(<f>2jx ~ Aj^2j - Wplj) = 0, 



which are satisfied if (3.1) holds. The remaining relation 

= [U,V( n) ]+Po 



-2<^ 2 ,0!> 
2<Vi,0 2 > 

is nothing but the constraint (3.2). Therefore the equation (3.4) holds if the system (3.1) holds together 
with the symmetry constraint (3.2). This means that and U constitutes a Lax pair for the constrained 
system consisting of (3.1) and (3.2). 

Lax operators can be represented explicitly. For example, in the case of n = 1, we have 

A(X) = -1 + Ef=i x^jjftMy " foifai) ■= -1 + A), 



B(A) = Ef=i 



Bn 



(3.5) 



C(A) = EjLi T^TT^i^ := Co. 



The Poisson bracket in this case is given by 



2 N 



(3.6) 



as defined in rcf. . For each case of n Jr Hffi use the standard Poisson bracket like the above. 

Now turn to constructing r-matrices&3 . A direct calculation, similar to refs.liS'G-il , gives rise to 

f {A(X), A(jj)} = {B(X), B(p)} = {C(X),C(fi)} = 0, 



- X-i 
_ Po 



Ml 



(B(Ai)-B(A)), 


(3.7) 


(C(A)-C( M )), 


(A(»)-A(X)), 





G 



for A(X),B(X),C(X) defined by (3.5). Setting 

M(A) = 



C(A) -A(X) 



, Mi = M (g> 7, M 2 = 7 ® M, P 



^ 1 \_ 

f 

f 

\ o o o f / 



we have 



{Mi(A),M 2 (/i)} 



/ {A(A),AO)} {A(A),BO)} {B(A),AO)} {B(A),BO)} \ 

{A(A),C(/i)} {A(A),-A(/i)} {B(A),CM} {5(A), 

{C(A),A( M )} {C(A), !?(/*)} {-A(A),A(m)} {-A(A),B(a*)} 

V {C(a),c( m )} {C(a),-4( m )} {-^(a),-^)} {-a(a),-a( m )} y 



Mi(A)+M 2 (/i) = 



[P,Mi(A)+M 2 (/z)] 



/A(A)+A(/i) B(/x) B(X) \ 

C(/i) A(A)-AO) B(X) 

C(A) -A(A)+A(/i) B(fi) 

C(A) C(/x) _A(A)-A(/x)y 

B(/z) - B(A) B(A) - BO) 

C(A)-CO) 2(-A(A)+AO)) BO)-B(A) 

CO) - C(A) 2(A(A) - A(/i)) B(A) - BO) 

C(A) - CO) CO) - C(A) 



Comparing (3.9) with (3. if) and using (3.7), we obtain the classical Poisson structuretH 

{Mi(A),M 2 0)} = [B(A,/i),Mi(A) + M 2 0)] 

with the following r-matrix 



B(A,/x) 



X — [i 



-P. 



It follows from (3.12) that 0-0 

{M| (A), M|0)} = [Bi(A ;M ),M 1 (A) +M 2 0)], 

= E<=oEj=o^W 1 - < M 2 (/i) 1 -^(A,/i)A^(A)A4(/i). 

This equality implies the relation 

{TrM 2 (A),TrM 2 0)} = 0. 

Therefore the integrals of motion generated from TrM 2 are in involution in pairs. 
Let us now show more examples. For n = 2, we have 



For n = 3, we have 



A(X) = -X + A Q , B(A) = 4>n+i + Bq, C{X)=^j n+1 +C . 



A(X) = -A 2 + \4> N+1 <t> N+2 + Aq, 



B(A) = (/) N+ iX + 2if> N+2 + B , 
C(A) = 0JV+2A - 2<ip N+1 + C . 
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For n = 4, we have 

A(X) = -A 3 + t}4>N+1^N+-A + ^<!>N+l1pN+2 + ^N+3<pN+2 + A), 

B(X) = (f> N +lX 2 + 07V+2A + <f> N +3 - + B , (3.17) 

t C(A) = 0at +3 A 2 + -0at +2 A + %jj N+ i - + Co- 

in the case of n = 3, the commutator relation (3.7) is established only for f3o — 1. In the cases of n = 2 
and n = 4, the commutator relation (3.7) is established only for /3o = — ^. All above result for n = 1 can 
be extended to the cases of n = 2, 3, 4. Therefore all corresponding constrained systems have the -(--matrices 
determined by (3.13). 

4 Integrals of motion and Liouville integrability 

First of all, we point out that if <j>ij,ipij, i — 1,2 (1 < j < N), are solutions of (3.1) with A = Xj, then a 
direct calculation yields ©H that 

^ Fi= 0, F,'= lM« + ihjfaj- (4- 1 ) 

From the relation 

which may be resulted directly from the Lax representation Nx = [U, N^], we can get another set of 
integrals of motion {F TO }£° for each case of n. Now we list them as follows. 

i) n = 1. The corresponding constraint is v =< ip2,<f>i >, w —< tpi,4>2 >> and the integrals of motion 
F„ has been found in refJ 10 ' 

Fx =< 1p 2 , 02 > - < ipi,<pi >, 

Fn z^i=l {^i^n—i "t - biCn—i) 2(X n: Ti ^ 2, 

where 

a; = i(< A l - X i\> x ,<$>x > - < A 1 - 1 ^,^ >), =< A^faAi >, q =< A^Vi,^ > ■ 

By using the idea in Refs.]l3t the required 2N integrals of motion for Liouville integrability of the 

constrained system in the case of n = 1 can be chosen as 

Fx,F2,-,Fn,Fi,-,F n . (4.2) 

ii) n — 2. The corresponding constraint is =< ip2,<j>i >, = — < ipi,4>2 >, and similarly the 
required 2iV + 1 integrals of motion can be chosen as 

Fi, Fa, Fn,F 2 , F 3 , .... Fv+2, (4.3) 
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where 

F 2 = y(<ipi,</>i> - <1p2,<h>) + <l>N+l1pN+l, 

f 3 = i(<,4Vi,'Ai>-<^ 2 ,^ 2 >)-^ti<Vi,^2>- :5 ^ ±1 <V2,0i>, h = -\f s , 

F n = \{< A n ~ 2 ipi,4>i >-< A n -^ 2l 4>2 >) - < ^"- 3 ^i,02 > -^tl < A»- 3 ^,0i > 

+ (aia k + bic k ), n > 4, 

;+fe=n,fe,;>2 

a ; = -|(< A'- 2 V>i,0i > - < A'- 2 ^ 2 ,^2 >), 6; = -5 < A l - 2 ^ 2l 4>i >, ci = -\< A l ~ 2 ipi, 2 > . 

iii) n = 3. The corresponding constraint is ^(v xx — 2u 2 w) =< ip2,<f>i >i \{ w xx — 2vw 2 ) =< tpi,(f>2 >■ 
The required 2N + 2 integrals of motion can be chosen as 

Fi, F2, Fat, F 3 , F 4 , -, Fv +4 , (4.4) 

where 

F 3 = - < Ipl, 01 > + < tp 2 , <j>2 > +2l/>Ar+20W+2 - 2tp N+1 (j) N+1 , 

F A = - < Alp!, fa > + < Alp 2 ,<h > + 0N+1 < ^1,02 > +4>N+2 < ^2, fa > +\<t>N +l<t>N +2 ~ 4lp N+1 1p N+2 = H, 
F 5 = - < A 2 1pi,fa > + < A 2 7/> 2 ,0 2 > + 0N+1 < Alp!, fa > + 0N+2 < Alp 2 ,fa > 

+ 2lp N+2 <l/j 1 ,<p2> -2^W+1 < "02, 01 > +\(f>N+l(t>N+2{< 1pl,<t>l> ~ < 1p2, 02 >), 
F n = -< A n ~ 3 fa, 01 > + < ^"- 3 ^2, 02 > +0JV+1 < A""Vl, 02 > +0JV+2 < A n ~^2, 01 > 

+ 2^ N+2 < A"" 5 Vl,02 > -2V>iV+l < ^"-^2,01 > +50W+10iV+2(< A"" 5 Vi,0i > 

- < A"- 5 ?A 2 ,0 2 >) + ^ (aia k + bic k ), n>6, 

2+fc=n,fc,Z>3 

a; = \{< A 1 - 3 ^, 0i > - < A'- 3 ?A 2 , 2 >), b x =< A l ~ 3 ip 2 , 0i >, q =< A'" 3 ^!, 2 > • 

iii) n = 4. The corresponding constraint is |(w IZI — 6wwi; x ) =< fa, fa >> j(wui — 6t)WMij;) = — < 
fa, fa > • The required 2 AT + 3 integrals of motion may be chosen as 

Fi, F 2 , Fjv, F 4 , F 5 , F 6 , Fv +6 , (4.5) 
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where 

Fa = \{< V'l, 01 > - < 1p2, 02 >) + \<t>N + 10W+3 + 4>N+l1pN + l + 0iV+3'0iV+3 + 4>N+2'4>N+2, 

F 5 = \{< AlPufa >-< AlP 2 ,fa >) - ^1 < lh,<h > -^T 1 < ^,01 > +0^+2^+1 
+ <t>N+3tpN+2 + %(f)N+l<f>N+2'lp% + 3 + %<p 2 N+1 1pN+2lpN+3, H = ~^F 5 , 

F 6 = J(< > - < ^V>2,02 >) - fa^+3 (< ^,0! > - < V>2,02 >) ~ ^ < ^/>1,02 > 

_3^±3 < Alht<h > _<h+l < ^ > -t^tl < ^ > +l(^ +1 ^ +2 + ^ +2 ^ + 3 
+20^+1*0^+20^+30^+2) + 0JV+3-0Af + l - g0J\r+30JV+l^Ar+3 ~ g0w+10JV+10J\T+3 + ^0^+10^+3, 

F 7 = \{< A 3 l[)i, 4>\ > - < A 3 ?/> 2 ,02 >) ~ \(t>N+l^N+z{< >ty>i,0i > - < A^ 2 ,02 >) 

-|(0W+10W+1 + 0^+3^+2 )(< 01>1 > - < 02,02 >) - < ^01,02 > 

_0JV±3 < ^ 0i > _ 0J^2 < HiA2 > _^N+1 < ^ 01 > 

-^(0w+3 - g0|r+l^iV+3) <ipi,<h> -\{i>N+i - §0^+1^+3) < 02, 01 >, 
= £(< ^"-Vl, 01 > - < A""V 2 , 02 >) - |0Ar+lV^+3(< A n "Vl, 01 > - < ^"^2, 02 >) 

-|(0JV+1VJV+1 + VjV+30iV+2)(< A n ~ 7 1pl,(f>l >-< A n - 7 ^ 2 ,4>2 >) - < A"" 5 ^!, 02 > 

_0JV±i < ^-5^ 2j01 > _0^+2 < A n-6^ h > _0JV+2 < A n-6^ fa > 

-5(0^+3 - g0ir+i0iv+3) < A n - 7 ^i,<^2 > -5(^+1 - ^jv+iV'Ir+s) < ^4"~ 7 02,0i > 
+ ^ (aia k + bic k ), n > 8, 

l+k=n,k,l>4 

a; = -g(< A'-Vi, 01 > - < ^ _4 02, 02 >), k = -\ < A'- 4 V>2, 01 >, d = ~\< ^"Vi, 2 > . 

The involution property {Fi,Fj} — {Fi,Fj} = can be easily checked and the involution property 
{Fi,Fj} = is a direct conclusion from the r-matrices in §3. Now the remaining problem for Liouvillc 
integrability is to prove that all the required integrals of motion are functionally independent. This can be 
achieved, indeed. The following analysis aims to providing a mathematical proof for it. 

i) The case of n — 1. Assume that the result on the functional independence is not true. Then there 
exist 2N integrals a\, ...a^,l3\, ...{3n, satisfying J2iLi a l + SjLi 0j ^ 0, such that 

JV JV 

i=i j=i 
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Since TUT = &3 S jl> 



dF n[ 



fe=o — -A™ tpu, 



( dFi 



A := dct 



dFi 



dct 



/ V>n 


V>21 

V o 



IplN 





27V 



dF n _ 

Q<foi 


A™ 1 V , 2;, n > 1, we have 


d Fn 

d<h 




dF N \ 

t)<h. 




d Fn 

d<j) 2 




dF N 
d<t> 2 I 


01=02=0 


—tbi i 

T -i- J- 


■■■ -Af-Vn N 




-V>1JV 
"021 


••• -Ajv 01AT 

••• Af"V 2 i 




4>2N 




1 4>2N J 




Ai • 


•• Af- 1 






A^v • 


\N-1 
■ ■ *N 







This means that all cti,/3j must be zero. Therefore the functions Fk,Fk(l < k < N) can be functionally 
independent at least on certain region of R 2N . 

ii) The case of n = 2. We need to show that the required 2N + 1 integrals of motion are functionally 
independent. Since we have 



8F 2 1 , dF 2 I 



dF 2 



dF n 

^jjl01=02=0jv + l = 

dF n 

^l=02=0N + l=O 



\ n— 3 



1 \n— 3„ 



30 



<?2i 

"T 2 — |01=02 = 0N + 1=O = 

JV+1 



4 A r 2 V>2;, n > 3, 



<9F, 



30 



IL -| ( ^ 1 =0 2 =0 N+1= o = 0, n > 3, 



N+l 
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we can make the following calculation 



A := det 



I dFi 

dFi 
dFi 



\ d(pN+i 
( o 



det 





^21 



V o 



IplN 





1p2N 





= (-1)^+1 



ii 



o 

^21 



= v>jv+i rifci v>2j 



<3<h 

OFn 
d<t>N+\ 



dF 2 
Ml 
dF 2 

pF 2 

d<t>N+i 



2^ 11 ' a: 



-^21 



V'W+l 





-^AiV^i 







-^Anfei 



arc 



W+2 



dF N+2 

C></>2 



jV+2 



»AT+1 



»1=¥>2=<?>JV + 1 = 



J - 2 AfV2i 



2 A 1 rll.x 



- 2 AfV2i 



1 \Nj, 



AfVn + |Af-V 1M 



1 \JV 

2 



^Aj^VliV + ^A^ l 1plN,x 



1 Ai 



AT-1 



1 A 



JV 



A 



JV-l 
N 



#0. 



Therefore the required integrals of motion are functionally independent at least on one region. 

iii) The case of n = 3. The required 2N + 2 integrals of motion are functionally independent, because we 
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have 



A := 



det 



( dFi 

Ufa 

dFi 

dFi 

dFi 
\ d(j) N+ 2 

( V'n 





$21 




V o 



Ipll 



$21 



OFn 
OFn 



OFn 
<> Fx 

d<t>N+2 d<pN+2 



dF 3 

~3fa 

dF 3 

m_ 

c>$v+i 

9^3 



dF N+4 \ 
~Ufa~ 

dF N+4 

8Fn+4 

CItpN+l 

<9F/v+4 

<3<t>N+2 ' 



4>i — 02 — o 

<t>N + l = <pN + 2 — i'N+l = 



-ipn -Ai^i 



-AfVn 







-V>1JV 
$21 



Ai^21 A?V21 + "2ip N+2 ip 
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Af +1 V2i + 2Af - 



11 



$2JV 1p2N Ajv$2iV X%1p 2 N + 2$v+2'$jV 



-4?/>7V+2 

2^+2 
-AfVn 

i>lN -X%1plN 

A?$n + 2ipN+2i>u 

lp2N >? N ^2N + 2lp N+2 ^iN 



A^ +1 $2W + 2\%- 1 ip N+2 A N 





-Af+Vi 



2 W+3 ^+2 nf=l V>li(A^2i + $V+2^) 



1 Ai 



-A^+Viiv 
Af +1 V2i + 2Af-Vjv+2^uv 

A^ +1 $2iV + 2A^ _1 $W+2$LiV 
A l 



1 A 



AT 



.JV-1 



^0. 



iv) The case of n = 4. The required 2iV + 3 integrals of motion are functionally independent, because we 
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have 



det 



( dFi 


dFw 


dF A 


dFN+e \ 


d<k ■ 








dFi 


d Fn 


dFi 


dF N +6 




d<p 2 


~d~$2~ 


d(/)2 


dFi 


dF N 


dF 4 


dFN+6 




dipN+i 


dipN+i 


dtpN+i 


dFi 


dF N 


dF 4 


dFN+6 


dtpN+2 




dlpN+2 


d4>N+2 


1 d Ft 


dF N 


dF 4 


dFN+6 




VlpN+3 


dlpN+3 


di>N+z ' 



det 



( ^11 


"021 





V o 



ryVll 





lp2N 







^lN 



■\^2N 




t>N+3 



-^1021 





<j>N+3 





-\\\^21 



— ^X%1p2N 
4>N+3 




1 \ N + 2 ,I, 



01—02 — 

0N + 1 — 0JV + 2 = 

IpN+l = V'IV + 2 = i>N + 3 = 



2^1^21 - 2^+3^11 



-g^AT^JV - ^N+^IN 







1 \ N + 2 nh 

■^^N YIN 
4< +2 ^21 - 



1 \iV+2 , 1 \N-li i 

~2^N V2N — 2 A iV (PN+3W1N 







1 Ax ••• Af 



iV-1 



1 A 



iV 



,iV-l 



#0. 



The above analysis allows us to conclude that the Hamiltonian systems (2.4), (2.9), (2.16) and (2.19) are 
completely integrable in the Liouville sensed . 



5 Conclusion and remarks 

We have successfully extended binary nonlinearization of AKNS spectral problem to higher-order symmetry 
constraints. For the first four orders n = 1,2,3,4 of the symmetry constraints, we explicitly presented 
the Hamiltonian structures, the Lax representations, r-matrices and independent integrals of motion in 
involution. Binary nonlinearization under the symmetry constraint in the case of n = 3 has also been 
considered by Xul£y but no Lax representations and r-matrix has been presented till now. 
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Indeed, for the case of n = 2k + 1, binary nonlinearization of AKNS spectral problem can always be 
successfully handled. However when n is an even integer greater than 4, we do not yet know any concrete 
expression for new dependent variables which put the constrained systems to be Hamiltonian. We hope that 
we can find a rule to guide us to introduce new dependent variables such that the constrained systems can 
be transformed into Hamiltonian forms. 

Throughout the paper, we suppose that the potentials v, w and their derivatives with respect to x tend to 
zero when \x\ — > oo. In the case without this zero boundary condition, the same problem is very interesting 
to be considered. Binary nonlinearization of 3 x 3 matrix spectral problemH' or even n x n matrix spectral 
problemEII for AKNS hierarchy under higher-order symmetry constraints is another problem which needs 
to be investigated. We will deliver some analysis to the problems in a separate paper. 

Acknowledgments 

One of the authors (Y. S. Li) is grateful to City University of Hong Kong for kind invitation and warm 
hospitality. This work is supported by the City University of Hong Kong, the National Basic Research 
Project of Nonlinear Science, and the Ministry of Education of China. 

References 



[i] 

[2] 
[3] 

[4] 

[5 

[G] 

m 

[8 

[9 
[10 

[11 
[12 

[13 
[14 



C.W.Cao: Chin. Q. J. Math. 3(1988) 90; Sci. China A (1990) 528 

C.W.Cao, X.G.Geng: J. Phys. A Math. Gen. 23(1990) 4117; J. Math. Phys. 32(1991) 2323 

Y.B.Zeng, Y.S.Li: J. Math. Phys 30(1989) 1671; Y.B.Zeng: Phys. Lett. A. 160(1991) 541; Y.B.Zeng, 
Y.S.Li: Chinese Adv. in Math. 24(1995) 111 

W.X.Ma: Nonlinear Physics, eds. C.H.Gu, Y.S.Li, G.Z.Tu (Springer Berlin 1990) p79; Acta Math. Appl. 
Sin. 9(1993) 92 

M.Antonowicz, SWojciechowski: Phys. Lett. A 147(1990) 455; J. Math. Phys. 33(1992) 211; Inverse 
Problem 9(1993) 201 

X.G.Geng, W.X.Ma: IL Nuovo Cimento A 108(1995) 477 

B.Konopelchenko, J.Sidorcnko, W.Strampp: Phys. Lett. A 157(1991) 17; Y.Cheng, Y.S.Li: Phys. Lett. 
A 157(1991) 22 

B.Konopelchenko, W.Srampp: Inverse Problem 7(1991) L17; J. Math. Phys. 33(1992) 3676; J.Sidorenko, 
W.Strampp: Inverse Problem 7(1991) L37 

B.Xu, Y.S.Li: J. Phys. A Math. Gen. 25(1992) 2957 

W.Oevel, W.Strampp: Commun. Math. Phys. 157(1993) 51 

W.X.Ma, W.Strampp: Phys. Lett. A 185(1994) 277 

W. X.Ma: J. Phys. So c. Japan 64(1995) 1085; Physica A 219(1995) 467; Chin. Ann. of Math. 18B(1997) 



79 ( Mv-mt/9512002| ) 



W.X.Ma, Q.Ding, W.G.Zhang, B.Q.Lu: II Nuovo Cimento B 111(1996) 1135 
W.X.Ma, B.Fuchsteiner, W.Oevel: Physica A 233(1996) 331 



15 



[15] L.D.Faddeev, L.A.Takhtajan: Hamiltonian Methods in the Theory of Solitons (Springer Berlin 1987) 

[16] Y.B.Zcng, J.Hictarinta: J. Phys. A Math. Gen. 29(1996) 5241 

[17] Y.B.Zcng: Inverse Problems 12(1996) 797; J. Math. Phys. 38(1997) 321 

[18] O.Babelon, C.M.Viallet: Phys. Lett. B 237(1990) 411 

[19] V.I. Arnold: Mathematical Methods in Classical Mechanics (Springer New York 1980) 
[20] X.X.Xu: J of Shandong Mining Institute 16(1997) 445 

[21] X.X.Xu: An nxn matrix spectral problem for AKNS hierarchy and its binary nonlinearization, preprint, 
1997 



16 



